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ABSTRACT 

From a quantitative analysis of 413 Virgo cluster early-type dwarf galaxies (dEs) with Sloan Digital 
Sky Survey imaging data, we find that the dE class can be divided into multiple subpopulations that 
differ significantly in their morphology and clustering properties. Three dE subclasses are shaped like 
thick disks and show no central clustering: (1) dEs with disk features like spiral arms or bars, (2) 
dEs with central star formation, and (3) ordinary, bright dEs that have no or only a weak nucleus. 
These populations probably formed from infalling progenitor galaxies. In contrast, ordinary nucleated 
dEs follow the picture of classical dwarf elliptical galaxies in that they are spheroidal objects and are 
centrally clustered like E and SO galaxies, indicating that they have resided in the cluster since a long 
time, or were formed along with it. These results define a morphology-density relation within the dE 
class. We find that the difference in the clustering properties of nucleated dEs and dEs with no or only 
a weak nucleus is not caused by selection biases, as opposed to previously reported suggestions. The 
correlation between surface brightness and observed axial ratio favors oblate shapes for all subclasses, 
but our derivation of intrinsic axial ratios indicates the presence of at least some triaxiality. We 
discuss possible interrelations and formation mechanisms (ram-pressure stripping, tidally induced 
star formation, harassment) of these dE subpopulations. 

Subject headings: galaxies: dwarf — galaxies: structure — galaxies: evolution — galaxies: statistics 
— galaxies: fundamental parameters — galaxies: clusters: individual (Virgo) 



1. INTRODUCTION 

As the most numerous type of galaxy in clusters, 
early-type dwarf galaxies are ideal probes to study the 
physical processes that govern galaxy formation and 
evolution in environments of different densit y. The 
pronounced morphology-de nsity relation (e.g.. iDressleil 
119801 : iBinggeli et al.l I1987D suggests that early-type 
dwarfs were either formed mainly in high-density envi- 
ronments, or originate from galaxies that fell into a clus- 
ter and were morphologically transformed. However, the 
actual f ormation mechanisms a re still a matter of de- 
bate (see lJerien fc Binggelill2005l and references therein). 
Most of the proposed scenarios are based on the vigor- 
ous forces acting w ithin a cluster environment, like ram- 
pressure stripping fCunn fc Gott''T972') of dwarf irregu- 
lar (dirr) galaxies (e.g.. van Z ee et al. 20041. tidally in - 
duced star formation in dirrs (fPavies fc Phillippg|ll988f ). 
or so-called harassment (jMoore et al.lll996l) of infalling 
late-type spirals through close encounters with massive 
cluster members. 

Early-type dwarfs form a rather heterogeneous class 
of o bjects. In addition to the classical dwarf ellipti- 
cals, ISandage fc BinggelH (|1984l ) introduced the class of 
dwarf SO (dSO) galaxies, which were conjectured to have 
disk components, based on signatures like high flattening 
or a bulge-f disk-like profile (Binggeli fc Cameron 1991). 
The identification of spiral substructure then provided 
the first direct pro of for a disk in an early-type dwarf 
(jJerien et al.ll2000l ). which, however, had not been clas- 
sified as dSO but as dwarf elliptical. Inspired b y simi- 
lar d i scoveries (e.g.. IBara zza et al. 2002; Graha m et al.l 
120031: iDe Riicke et al.l 12*0031) , we performed a search for 



disk features in 41 Virgo cluster early-type dwarfs 
(Lisker et al."2006a, hereafter Paper I). We thereby in- 
cluded galaxies classified as dwarf elliptical and as dSO 
to avoid any preselection bias, and assigned them the 
common abbreviation "dE" , which we adopt for this Pa- 
per as well. We identified disk features in 36 dEs, and 
argued that they constitute an unrelaxed population of 
disk-shaped galaxies different from the classical dwarf el- 
lipticals (Paper I). 

But the dE class shows yet more diversity: nucleated 
and n on-nucleated dEs ha.ve different clustering prop - 
erties (van den Berghlll986l : iFcrgvison fc Sandage"l989t). 
their fi attening distribution s differ ( Binggeli fc Cameron 
1991t iRvden fc Terndrupl [19941 : iBinggeh fc Popescu 
1995D . and color differences were reported as well 



lijjjoil. ana color a iiteren ces were rep ortea as wei 
(|Rakos fc Schomberlj [200l iLisker et"alf I2005D . More- 
over, several of the bright dEs display blue central regions 
caused by recent or ongoing star formation (Lisker ct aQ 
l2006bl hereafter Paper II), and also differ in their spatial 
and flattening distributions from the bulk of dEs. Thus, 
prior to discussing possible formation mechanisms, we 
need to systematically disentangle the various dE sub- 
classes observationally. This is the purpose of this Paper. 

2. SAMPLE SELECTION 

While our dE sample selection was already described 
in Papers I and II of this series, these studies were still 
based on the Data Release 4 of the Sloan Dig ital Sky 
Survev (SDSS: lAdelman-McCarthv et"al]l200l . Since 
we are now using the full SDS S Data Release 5 (DR5, 
lAdelman-McCarthv et al.ll2007| ) dataset, we provide here 
a detailed, updated description of our selection. 
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Th e Virgo Cluster Catalog (VCC, iBinggeli et all 
198511 ■ along with revised classifications from 
Barazza et air(l2002l VC C 142 2). iBarazza et alj ()2003l 
VCC 0850), and iGeha et"aI1 (|2003|, VCC 1488), con- 
tains 1197 galaxies classified "dE" or "dSO", including 
candidates, that are certain or possib le cluster mem- 
bers a ccording to IBinggeli et al.l (|1985f) . IBinggeli et al.l 
()1993l ). and Paper II. 552 of these fall within our 
chosen limit in apparent B magnitude from the VCC 
of niB < 18.0 mag (see Paper I). This is the same 
magnitude li mit up to which the VCC was found to 
be complete (|Binggeli et al.l Il985| ) . When adopting a 
Virgo cluster distance oi d = 15.85Mpc, i.e., a distance 
modu lus m — M = 31.0 mag (see, e.g., iFerrarese et al.l 
l2000l) . which we use throughout, this corresponds 
roughly to a limit in absolute magnitude of < —13.0 
mag. 

Six galaxies are not covered by the SDSS. While we ini- 
tially included objects with uncertain classification (e.g., 
"dE?"), we then excluded all 50 galaxies that appeared 
to be possible dwarf irregulars from visual inspection of 
the coadded SDSS g,r, and i images (see Paper I), or were 
classified as "dE/Im". Three more objects (VCC 0184, 
0211, and 1941) were excluded because they appear to 
be probable background spirals. Finally, VCC 1667 could 
not be classified properly, since it is significantly blended 
with multiple other galaxies. This leads to a final dE 
sample of 492 certain or possible cluster members, con- 
taining 426 certain cluster members on which we focus 
in the present Paper. 

2.2. Presence of nuclei 

While our classification of nucleated and non-nucleated 
dEs relies on the VCC, it is known from HST observa- 
tions that many apparently non-nucleated dEs actually 
host a fai nt nucleus hardly detecta ble with ground-ba sed 
imaging (|C6te et al.ll2006L also see lLotz et al.ll2004b| ) . A 
direct cor nparison of the VC C classification with the re- 
sults from lCote et al.l ()2006D shows that, as a rough rule 
of thumb, the detection of dE nuclei in the VCC be- 
comes incomplete for nucleus magnitudes that are fainter 
than the respective value of the host galaxy's central 
surface brightness, measured within a radius r = 0.1" 
(Fig.[T]). Our non-nucleated dEs could thus be more ap- 
propriately termed dEs without a nucleus of significant 
relative brightness as compare d to the underlying light 
of the galaxy's center. In fact, iGrant et al.l ()2005D sug- 
gested that dEs classified as nucleated and non-nucleated 
might actually form a continuum of dEs with respect to 
relative nucleus brightness. Therefore, the VCC classi- 
fication basically translates into probing opposite sides 
of this continuum — and this is exactly what makes it 
useful for our study of dE subclasses. If the relative 
brightness of a nucleus depends on its host galaxy's evo- 
lutionary history, then one might expect nucleated and 
"VCC-non- nucleated" dEs to exhibit different popula- 
tion properties. 

3. DATA 

The SDSS DR5 covers all VCC galaxies except for an 
approximately 2° x 2?5 area at a « 186?2, i5 ^ +5°0. 
It provides reduced images taken in the u, g, r, i, and 
z bands with an effective exposure time of 54s in each 
band (see also lStoughton et "al]|2002f ) . as well as the nec- 



essary parameters to flux calibrate them. The pixel scale 
of 0'.'396 corresponds to a physical size of 30 pc at our 
adopted Virgo cluster distance of d = 1 5.85Mpc. The 
SDSS imaging camera ()Gunn et al .11 19981 ) takes data in 
drift-scanning mode nearly simultaneously in the five 
photometric bands, and thus combines very homoge- 
neous multicolor photometry with large area coverage 
and sufficient depth to enable a systematic analysis of 
dEs. The images have an absolute astrometric accuracy 
of rms < (f.'l per coordinate, and a relative accuracy 
between the r b and and each of the other bands of less 
than 0.1 pixels (jPier et al.ll200l . They can thus easily 
be aligned using their astrometric calibration and need 
not be registered manually. 

The rms of the noise per pixel corresponds to a surface 
brightness of approximately 24.2 mag arcsec"^ in the u- 
band, 24.7 in g, 24.4 in r, 23.9 in i, and 22.4 in z. The 
typical total signal-to- noise ratio (S/N) of a bright dE 
i'm-B.vcc ~ 14) amounts to about 1000 in the r-band 
within an aperture radius of approximately two half-light 
radii. For a faint dE (me ~ 18) this value is typically 
about 50. While the S/N in the g and i-band is similar, 
it is several times lower in the z-band and more than ten 
times lower in the u-band. 

The SDSS provides photometric measurements for our 
galaxies, but we fou nd these to be incorrect in many cases 
(jLisker et al.ll2005h . The SDSS photometric pipeline sig- 
nificantly overestimates the local sky flux around the 
Virgo dEs due to their large apparent sizes and low sur- 
face brightness outskirts. This affects the derivation of 
isophotal and Petrosian radii, the proflle flts, and sub- 
sequently the calculation of total magnitudes, which can 
be wrong by up to 0.5 mag. For this reason, we used B 
magnitudes from the VCC throughout the first two pa- 
pers of this series. In the meantime, we have performed 
our own structural and photometric measurements (see 
Sect. 2]), which we shall use here as well as in future pa- 
pers of this series. Still, when we refer to B magnitudes, 
these were adopted from the VCC. 

Heliocentric velocities for part of the sample are 
provided by the NASA/IPAC Extragalactic Database 
(NED; also see Paper II for more detailed references). 

4. IMAGE PREPARATION AND ANALYSIS 

4.1. Sky subtraction 

The sky level on the SDSS images can vary by some 
tenths of the noise level across an image. For a proper 
determination of Petrosian radii of the dEs (see Sect. 14. 3( ) 
despite their low surface brightness outskirts, it is thus 
not always sufficient to subtract only a single sky fiux 
value from each SDSS image. Therefore, we performed 
sky subtraction through the following procedure. First, 
we constructed object masks for each SDSS image from 
the so-called segmentation images of the Source Extrac- 
tor software (|Bertin fc Arnoutil 19961 ) by expanding these 
through smoothing with a Gaussian filter (using IRAF^, 
l^dy 1993). A preliminary sky level was then determined 
for each image as the median of all unmasked pixels, 
clipped three times iteratively at 3cr. In order to reach a 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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higher S/N than that of the individual images, we then 
produced a coadded image by summing the (weighted) 
g, r, and i-band images as described in Paper I. We then 
obtained an improved object mask from the coadded im- 
age and used this to refine our sky level measurement. 

Finally, the sky flux distribution across the image was 
determined by computing the average flux - clipped 
flve times iteratively at Str - of all unmasked pixels in 
201x201 pixel boxes, centered every 40 pixels. This grid 
of values can be stored as a 52x38 pixel "sky image". 
Pixels in this sky image that did not contain useful val- 
ues due to too many masked pixels in the parent image 
were linearly interpolated using IRAF fixpix. We then 
applied a 3x3 pixel median filter to the sky image, ex- 
panded it to match its parent SDSS image's size (using 
IRAF magnify with linear interpolation) , and subtracted 
it from the latter. This yields the final u, g, r, i, and z 
images. 

We point out that there is, to our knowledge, no gen- 
eral agreement or recipe as to whether to use, e.g., the 
clipped mean, the median, the clipped median, or the 
mode, for determination of the sky level. However, it is 
advisable that the chosen approach be reconciled with 
the image measurements to be performed, which in our 
case is the derivation of Petrosian radii (see Sect. 14. 3[) . 
Since the latter is based on the average flux within given 
annuli, we chose to use the clipped average flux of all 
unmasked pixels for our sky level measurement. This 
guarantees that the resulting flux level in each image is 
zero as "seen" by the Petrosian radius calculation.^ 

4.2. Calibration and extraction 

We calibrated the sky subtracted SDSS images using 
the provided flux calibration information (photometric 
zeropoint and airmass correction). We also corrected 
for the reported SDSS zeropoint offse ts in the u and 
z bands from the AB system (Oke fc GunnI Il983l 
see http: / /www. sdss.org/ dr5 / algorithms /fluxcal.html) . 
However, before working with the images, it is advan- 
tageous to put together adjacent images: a number of 
galaxies partly extend beyond the image edges and reap- 
pear on the corresponding neighbouring image. Bright 
dEs typically have apparent diameters of 300 pixels 
or more, which is rather large compared to the SDSS 
image size of 2048x1489 pixels. The SDSS astrometric 
calibration allows us to accurately put together adjacent 
images, which we did before extracting an 801x801 
pixel cutout image for each galaxy. These cutout images 
were then corrected for Galactic extinction, using one 
value per image, calcul ated with the dust rn aps and cor- 
responding software of 'Schlegel et al.' (l99§, provided at 
[http: / /www. astro. princeton.edu/^schlegel/dust/data/). 
From the g, r, and i cutout image we produced a final 
coadded image for each galaxy. 

4.3. Morphology 

^ The reason why such considerations are at all necessary is the 
same as that for which the SDSS pipeline overestimated the local 
sky flux: the Virgo dEs are large in apparent size and cover 10^ to 
10 pixels, but their low surface brightness outskirts cause a large 
number of these pixels to have S/N< 1. Thus, a wrong sky level 
estimate of the order of just a few tenths of the noise level can have 
a large effect in total. 



We perform an iterative process of determining shape 
and total flux for each galaxy, as described below. 
Throughout this process, we mask disturbing foreground 
or background objects, i.e., we do not consider masked 
pixels in any calcula tion. We star t with deriving 
the Petrosian r a dius fjPetrosianI Il976( l. as defined by 
IStoughton et al.l (|2002[ ). on the coadded image. Using a 
circular aperture with one Petrosian radius, we then find 
the center of the galaxy's image by itera tively searching 
for th e minimum asymmetry, following ' Conselice et al.l 
(|2000| ). The asymmetry A is calculated as 

i/i — /j,i8oi 



where fi is the flux value of the i-th pixel, and /i^iso is 
the flux value of the corresponding pixel in the 180-degree 
rotated image. 

The asymmetry is computed using an initially guessed 
central position (from Paper I for objects in the SDSS 
DR4, and from visual examination for obje cts in DR5, 
using SAOImage DS9. lJove &: Mandell[20M l. as weU as 
for using the surrounding eight positions in a 3 x 3 grid as 
center. If one of the surrounding positions yields a lower 
asymmetry, it is adopted as new central position. This 
process is repeated until convergence. We perform two 
of these "asymmetry centerings" : a flrst one with a step 
size of 1 pixel, and a second one with a step size of 0.3 
pixels. The initial and final value typically differ by less 
than a pixel. 

We then compute the parameters defining an elliptical 
aperture ( axial ratio and positi on angle) from the image 
moments ([Abraham et al.lll994| ) , and derive a "Petrosian 
semimajor axis" ap, i.e., we use ellipses instead of cir- 
cles in the calcul ation of the Petrosian radius (see, e.g., 
iLotz et al1l2004a[ ). Within this elliptical aperture with a 
semimajor axis a of lap, we perform another iteration to 
re-derive the elliptical shape parameters from the image 
moments, and also to re-derive ap. 

The elliptical shape is then applied to measure the total 
flux in the r band within an elliptical aperture with a — 
2ap, which also yields a half-light semimajor axis in r 
(ahi,r)- Using this value for ahi,r, we go back to the 
coadded image and flt an ellipse to the isophotal shape 
of the galaxy at a = 2ahi,r, using IRAF ellipse. The 
elliptical annulus used for the isophotal fit ranges from 
2°-''^au^r to 2i-25ahi,r. 

This new elliptical shape is now used to derive the final 
Petrosian semimajor axis on the coadded image, and to 
subsequently measure again the total flux in the r band 
within a = 2ap, yielding the flnal value for ahi,r- The 
isophotal shape is then measured again at a = 2ahi,r, 
yielding the axial ratio that we shall use throughout this 
Paper. 

Since we masked disturbing foreground or background 
objects by not considering their pixels, our measured to- 
tal flux for a given galaxy is always lower than it would be 
without any such "holes" in the galaxy's image. In order 
to correct for this effect, we subdivide the final aperture 
of each galaxy into 20 elliptical annuli of equal width, 
and assign each masked pixel the average flux value of 
its respective annulus. This yields our final value for the 
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total r band flux and the corresponding magnitude. The 
difference to the uncorrected value is typically less than 
0.1 mag. 

For 13 of our dEs, the derivation of the Petrosian radius 
did not converge, due to the fact that these galaxies sit 
within the light of nearby bright sources. While in some 
of these cases, it would still be possible to "manually" 
define an axial ratio for the galaxy, we decided to exclude 
these objects from our sample, since no reliable r band 
magnitudes can be derived, which are needed for our 
definitions of dE subclasses in Sect. [5] This leaves us 
with a working sample of 413 Virgo cluster dEs. 

5. EARLY-TYPE DWARF SUBCLASSES 

5.1. Subclass definitions 

Of our 413 Virgo dEs, 37 display disk features, like 
spiral arms, bars, or signs of an edge-on disk (Paper I, 
adding VCC 0751 to the objects listed there in order to 
update to SDSS DR5). We term these objects "dE(di)s", 
and separate this dE subclass from the ordinary, "fea- 
tureless" dEs (Fig. (2). In order to further explore the di- 
versity of the latter, we perform a secondary subdivision 
into nucleated ("dE(N)") and non-nucleated ("dE(nN)") 
galaxies, based on the identification of nuclei in the VCC 
as outlined in Sect. 12.21 Since a further subdivision of 
the dE(di)s would lead to statistically insignificant sub- 
samples, we shall instead discuss their nucleated fraction 
in the text. Finally, since our galaxies span a range of al- 
most 5 mag in r, it appears worth performing a tertiary 
subdivision into dEs brighter and fainter than the me- 
dian r brightness of our full sample, namely — 15.67 
mag. Moreover, all but three of the dE(di)s are brighter 
than this value; thus our subdivision allows us to com- 
pare them to ordinary dEs of similar luminosities. The 
percentage of each subsample among our full sample of 
413 dEs is given in parentheses in Fig. [2l whereas the 
actual number of galaxies contained in each subsample 
is given in the left column of Fig. [3l 

The subclasses defined so far are based on struc- 
tural properties only — for morphological classification 
of galaxies, it is not advisable to use color informa- 
tion. However, in Paper II we identified a significant 
number of dEs with blue centers (17 galaxies, including 
VCC 0901 from the SDSS DR5). These objects, termed 
"dE(bc)s", exhibit recent or ongoing central star forma- 
tion, similar to NGC 205 in the Local Group. They were 
morphologically classified as dwarf ellipticals or dSOs by 
ISandage fc Binggelil (|1984l ). and their regular, early- type 
morphology was confirmed in Paper II; thus, they are 
not possible irregular galaxies, which we have excluded 
from our samples here and in previous papers of this 
series. The fiattening distribution of the dE(bc)s was 
found to be incompatible with intrinsically spheroidal 
objects (Paper II), and their distribution with respect 
to local projected density suggests that they are an un- 
relaxed population. The latter result is similar to the 
spatial distribution of Virgo and Fornax dwarfs with 
early-type morphol ogy that are gas-rich and/or show star 
formation ( Drinkw ater et al.]l2001t IConselice et al.l[2003t 
iBuvle et al.i,2005.) . 

While it is not clear a priori that any of the dE sub- 
classes defined above are evolutionary interrelated, each 
dE(bc) unavoidably evolves into one of the above dE 
types once star formation ceases and the central color 



reddens (Paper II). Therefore, and because the dE(bc)s 
are defined through color instead of morphological prop- 
erties, we do not consider them a morphological dE sub- 
class.'^ On the other hand, their star formation and pres- 
ence of gas (Paper II) might imply that their formation 
process is not completely finished yet. It thus appears 
more cautious to separate them from the rest of dEs (see 
Fig. [2]) in order to not bias the population properties 
of the other subclasses. In the discussion (Sect. [7|) we 
try to assess which dE type(s) the dE(bc)s could pos- 
sibly evolve into. Note that four objects are common 
to both the dE(di) and the dE(bc) sample. We exclude 
these from the sample of dE(di)s, which now comprises 
33 galaxies. Table [1] lists our dEs along with their sub- 
class. 

A similar subdivision of the dE class into bright 
and faint (non-)nucleated s ubsamples was performed by 
[Ferguson fc Sandaed (|1989f ). also with the aim of study- 
ing shapes and spatial distributions of the resulting sub- 
sam ples. Our s ubdivision is different in two respects: 
first. [Ferguson fc Sandage defined all galaxies with tob < 
17.5 mag as "bright", whereas our magnitude separation 
(at rriy = 15.67 mag) is done at significantly brighter 
values and divides our full sample into equally sized 
halves. Second, we have the advantage of excluding 
dE(di)s and dE(bc)s from the "normal" dEs, thereby 
obtaining cleaner subsamples, especially for the bright 
objects: all but three of the dE(di)s are brighter than 
— 15. 67 mag. 

While iFerg uson fc Sandagd ()1989D found statistically 
significant differences in the spatial distributions of their 
subsamples - with dE(N)s being much more centrally 
clustered than the bright dE(nN)s - their flattening dis- 
tributions were only based on eye-estimated axial ratios 
from photographic plates. These can be uncertain by 
^20% (Ferguson fc Sandage 1989). With our measured 
axial ratios from the coadded SDSS images at hand, we 
therefore present in the following subsection a more de- 
tailed and accurate study of the flattening distributions 
of the different dE subsamples, and attempt to deduce 
their approximate intrinsic shapes. 

5.2. Subclass shapes 

From the axial ratio measurements of our galaxies 
(Sect. 14. 3p . we put together the flattening distributions 
of each dE subsample. These are presented in the sec- 
ond column of Fig. [3] as running histograms, i.e., at each 
sampling point we consider the number of objects within 
a bin of constant width, and normalize the resulting 
curve to an area of 1. The bin width is 0.15, which 
we have chosen to be one fifth of the range in axial ra- 
tio covered by our galaxies. The sampling step is 0.04 
(one quarter of the bin width). The bright and faint 
dE(N)s, and also the faint dE(nN)s, predominantly have 
rather round apparent shapes, while the bright dE(nN)s, 
dE(di)s, and dE(bc)s exhibit a significant fraction of ob- 
jects with rather flat apparent shapes. 

^ A morphological peculiarity of several dE(bc)s is that they 
show central irregularities, which are presumably due to gas, dust, 
and/or star formation, similar to NGC 205. These can be seen, e.g., 
when constructing unsharp mask images (Paper II). However, an 
attempt to quantify these weak features through image parameters 
like asymmetry or dumpiness yielded no clear separation from the 
bulk of dEs. Moreover, not all dE(bc)s display such features. 
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Since the division between bright and faint objects at 
TTir — 15.67 mag is somewhat arbitrary, we test whether 
the difference between the axial ratio distributions of 
faint and bright dE(nN)s becomes even more pronounced 
if a wider magnitude separation is adopted. The grey 
curves in the respective panels of the second column of 
Fig. [3] show the distributions for bright dE(nN)s with 
rrir < 15.67 — 0.5 mag (23 objects) and for faint dE(nN)s 
with rrir < 15.67 + 0.5 mag (86 objects). While the faint 
dE(nN)s basically remain unchanged, the bright dE(nN)s 
indeed tend slightly towards flatter shapes, but the dif- 
ference is rather small. 

A statistical comparison of the axial ratio distributions 
of our dE subsamples confirms what is seen in Fig. [31 a K- 
S test yields very low probabilities that any of the "flat- 
ter" subsamples (bright dE(nN)s, dE(di)s, and dE(bc)s, 
lower three rows) could stem from the same true dis- 
tribution function as any of the "rounder" subsamples 
(bright and faint dE(N)s as well as faint dE(nN)s, upper 
three rows). This confirms our findings from Papers I 
and II for the dE(di)s and dE(bc)s, respectively. The 
resulting probabilities from the K-S test for the pairwise 
comparison of the subsamples are given as percentages 
in Fig. m Interestingly, the lowest probability of all com- 
parisons is obtained when matching the distributions of 
bright and faint dE(nN)s: here, the probability of the 
null hypothesis that they stem from the same underlying 
distribution function is only 0.10%. Note that the prob- 
abilities for the comparison of the "flatter" subsamples 
with the "rounder" ones increase slightly with decreas- 
ing sample size, going from the bright dE(nN)s to the 
dE(di)s and then to the dE(bc)s. However, the probabil- 
ity for a common underlying distribution of dE(bc)s and 
the bright and faint dE(N)s is still only 3.8% and 4.6%, 
respectively. 

Is it possible to deduce the distributions of intrinsic 
axial ratios from t hose of the apparent ones? As dis- 
cussed in detail by iBinggeli fc Popescul (|1995f ). the in- 
trinsic shapes can be deduced when assuming that they 
are purely oblate or purely prolate. The distribution 
function 5" of intrinsic axial ratios q can then be derived 
from the di stribution function $ of observed axial ratios 
p through (jFall fc Frenldll983L eqs. (6) and (9)) 



*(9) = -v/r^;f fdp 
TT dq Jq 



$(p) 



for the oblate case, and 



*(?) 



2 



dq Jo 



(2) 



(3) 



for th e prolate case. Following Bing geH fc Popescul 
(|1995f ). we first defined adequate analytic functions $(p) 
that represent the observed distributions, and then eval- 
uated the above equations numerically. The analytic 
"model functions" are shown in the third column of 
Fig. [3l they were constructed from combinations of 
(skewed) Gaussians with each other and, in some cases, 
with straight lines. Note that, for the dE(bc)s, we de- 
cided not to follow the observed distribution in all detail, 
since it is drawn from a rather small sample of 17 galax- 
ies, which probably is the cause of the fluctuations seen. 

The deduced intrinsic distributions are presented in 
the fourth column of Fig. [3l for the oblate (grey lines) 



and prolate (black dash-dotted lines) case. We also show 
3-D illustrations of the galaxy shapes for each distribu- 
tion (fifth column), using in each case the axial ratio 
of the 25th percentile (left 3-D plot) and the 75th per- 
centile (right 3-D plot). These results confirm that the 
bright dE(nN)s, dE(di)s, and dE(bc)s do have lower ax- 
ial ratios than the bright and faint dE(N)s and the faint 
dE(nN)s. Furthermore, we point out that the bright and 
faint dE(N)s span a rather wide range of intrinsic axial 
ratios, and are, on average, somewhat flatt er than what 
was deduced bv iBinggeU fc Popescul (|1995D : our median 
value (see the 3-D illustrations in Fig. is slightly flat- 
ter than E3 for the prolate case, and slightly flatter than 
E4 for the oblate case. 

Can we decide whether the true shapes of our galax- 
ies are more likely to be oblate or to be prolate? For 
this purpose, we make use of the surface brightness test 
(iMarchant fc O lson'H^. 'R ichstonelll97l , again foUow- 
ing IBinggeli fc Popcscu (199^. If dEs were intrinsically 
oblate spheroids, galaxies that appear round would be 
seen face-on and should thus have a lower mean sur- 
face brightness than galaxies that appear flat; the latter 
would be seen edge-on. For the prolate case, the inverse 
relation should be observed. However, before we can per- 
form this test, we need to take into account the strong 
correlation of dE surface br ightness with magnitude (e.g., 
IBinggeli fc Cameronlll991l ): if, by chance, the few appar- 
ently round galaxies in one of our smaller subsamples 
would happen to be fainter on average than the appar- 
ently flat ones, this could introduce an artificial relation 
of axial ratio with surface brightness. Therefore, instead 
of directly using surface brightness like earlier studies 
did, we use the surface brightness offset from the mean 
relation of surface brightness and magnitude. We plot 
these values, measured in the r band within a = 2ahi,r, 
against axial ratio (measured at the same semimajor axis, 
see Sect. 14. 3p for each dE subsample, shown in the right- 
most column of Fig. [31 For all subsamples, a positive 
correlation of surface brightness offset with axial ratio 
can be seen, favoring the oblate model in agreemen t with 
earlier studies (e .g., Marchant fc Olson 1979; Rich ston"i 
ll979HBinggen fc Popesculll995l) . For the "rounder" sub- 
samples (top three rows), some additional contribution 
by prolate objects might be "hidden" within the rather 
large scatter of surface brightness offsets at larger ax- 
ial ratios. We denote these results in Fig. |3| by the ar- 
rows pointing from the surface brightness test diagram 
towards the favored intrinsic galaxy shapes. The arrow 
size represents the implied contribution from intrinsically 
prolate and oblate objects. Among the "flatter" sub- 
samples (lower three rows), for which the oblate case is 
favored, the dE(di)s have the lowest axial ratios, with 
a median value of 0.33 (bright dE(nN)s: 0.42, dE(bc)s: 
0.44). The galaxies in these subsamples are thus most 
likely shaped like thick disks. 

The above considerations needed to be restricted to 
purely oblate and purely prolate shapes. However, for all 
subsamples, a small part of the deduced (and favored) in- 
trinsic oblate distribution becomes negative at large axial 
ratios, trying to account for the low number of apparently 
round objects. This implies that most of the galaxies 
might actually have tri axial shapes, in accor dance with 
the conclusions of .Binggeli fc Popescul (|1995f) . 
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5.3. Subclass distribution within the cluster 

While it has been known for some time that nucleated 
and non-nuc leated dEs have diff e rent clustering prop- 
erties (e.g., Ivan den Berghl Il986t iFerguson fc Sandagg 

mthis statement has been challenged bv lCote et al. 
, who conjectured that it might just be the result 
of a selection bias in the VCC. It therefore appears worth 
to perform a quantitative comparison of the distributions 
of our dE subsamples within the cluste r, and to th e n pro - 
ceed with testing the issues raised by ICote et al] ()2006D 
in detail. 

The projected spatial distributions of our subsamples 
are shown in the middle column of Fig \5\ While both 
bright and faint dE(N)s exhibit a rather strong central 
clustering, the faint dE(nN)s appear to be only moder- 
ately clustered, and the dE(di)s and dE(bc)s show ba- 
sically no central clustering. The bright dE(nN)s even 
seem to be preferentially located in the outskirts of the 
cluster. 

To put the above on a more quantitative basis, we 
present in Fig. [S] the cumulative distribution of each of 
our subsa mples wi t h resp ect t o local projected den sity. 
Following 'Dressier (1980") and lBinggeli et al.l ()1987D . we 
define the latter for each galaxy as the number of objects 
per square degree within a circle that includes the ten 
nearest neighbours, independent of galaxy type. Only 
certain cluster members are considered. For comparison, 
we also show the same distributions for different Hubble 
types (Fig. [H inset), i.e., for the rather strongly cen- 
trally clustered giant early- type galaxies, as well as for 
the weakly clustered a nd probably infallin g spiral and 
irregular galaxies (e.g.. iBinggeh et al1ll987D . 

As a confirmation of the impression from the spa- 
tial distribution, the bright dE(nN)s are preferentially 
found in regions of moderate to lower density, similar 
to (and at even slightly lower densities than) the dis- 
tribution o f irregular galaxies, i n acco rdance with the 
findings of iFerguson fc Sandagd (|1989D . This imphes 
that they, as a population, are far from being virial- 
ized. The densities then increase slightly going from the 
bright dE(nN)s to the dE(bc)s, dE(di)s, and the faint 
dE(nN)s, in this order. Still, all of these are distributed 
similarly to the irregular and spiral galaxies in the clus- 
ter, again implying that they are unrelaxed or at least 
largely unrelaxed galaxy populations, and confirming the 
impression from their projected spatial distribution. In 
contrast, both bright and faint dE(N)s are located at 
larger densities, and display a distribution comparable 
to the E and SO galaxies, i n agre ement with the re- 
sults of IFerguson fc Sandagd 11989) . This would sug- 
gest that they are a largely relaxed or at least partially 
relaxed population. Note, however, that the Es alone 
(without the S O s) are located at still higher densities. 
IConselice et al.l ()2001[ ) pointed out that only the Es ap- 
pear to be a relaxed galaxy population, while all others, 
including the SOs, are not — thus, the dE(N)s presum- 
ably are not fully relaxed either. 

We performed statistical pairwise comparisons of the 
distributions of our dE subsamples with respect to den- 
sity, similar as for the axial ratios in Sect. 15.21 The K-S 
test probabilities for the null hypothesis that two ob- 
served distributions stem from the same underlying dis- 
tribution are given as percentages in Fig. [T] Even though 



the faint dE(nN)s are, among the "lower-density" sub- 
samples, closest to the bright and faint dE(N)s, their 
probability for having the same underlying distribution 
is 0.08 and 0.07%, respectively. These probabilities are 
higher for the dE(di)s and dE(bc)s: although they are 
located at even lower densities, their rather small sam- 
ple sizes let the probability increase as compared to that 
of the faint dE(nN)s. Finally, the bright dE(nN)s are 
located at such low densities that their K-S test compar- 
ison with the dE(N)s yields a probability of 0.00%, and 
that even the comparison with the faint dE(nN)s only 
yields a probability of 3.7% for them having the same 
true distribution. Given the morphological differences 
between the subsamples, as deduced in Sect. 15.21 Fig. [H] 
basically shows a morphology-density relation within the 
dE class. 

This view appears to be corroborated by the distribu- 
tions of heliocentric velocities (right column of Fig. [5]) 
of the dE subsamples: that of the bright dE(N)s has 
a single peak and is fairly symmetric, while especially 
the faint dE(nN)s, dE(di)s, and dE(bc)s display rather 
asymmetric distributions with multiple peaks. The lat- 
ter could be interpreted as being a signature of infallin g 
populations (jTullv fc Shavalll984Hc"onselice et al.ll2001h . 
However, the differences between these velocity distribu- 
tions are not or only marginally significant — the "most 
different" pair of distributions according to the K-S test 
are the bright dE(nN)s and the dE(bc)s, which have a 
probability of 6.6% for the null hypothesis. The main is- 
sue here are the small sample sizes: only a fraction of the 
galaxies of each subsample has measured velocities (num- 
bers are given in parentheses in the left column of Fig. [5]), 
which are available from the NED for 193 of our 413 dEs, 
and, e.g., for only 19 of our 39 bright dE(nN)s. Similarly, 
measurements of the skew or kurtosis of the distributions 
do not yield values that differ significantly from zero. We 
can thus only state that the rather asymmetric, multi- 
peaked distributions of the faint dE(nN)s, the dE(di)s, 
and the dE(bc)s would be consistent with our above con- 
clusion that they are mostly unrelaxed populations, but 
that more velocity data is needed to perform a reliable 
quantitative comparison of velocity distributions. 

5.4. Remarks on possible selection biases 

The different spatial distribution of dE(N)s and 
dE(nN)s was long considered a fundamental and well- 
fo unded obseryatioii , but has recently been questioned 
bv lCote et al.l ()2006f ). These authors argued that galax- 
ies with high central surface brightness (HSB, with 
Mg. central ^ 20 mag arcscc"^ or B < 14.55) would 
have been preferentially classified as non-nucleated in 
the VCC, which may have lead to a selection bias in 
the VCC that artificially relates spatial distribution to 
nucleus presence. We test this conjecture by considering 
the following points: 

(1) If the dE(nN)s were objects in which nuclei have 
preferentially gone undetected due to a too large cen- 
tral surface brightness, the dE(nN)s' surface brightnesses 
should, on average, be significantly higher than those of 
the dE(N)s. However, the mean surface brightness in r 
within the half-light aperture has very similar median 
values for the bright dE(nN)s {fir — 22.65 mag arcsec^^) 
and the bright dE(N)s {fir — 22.63 mag arcsec"^), which 
makes such a bias unlikely. Furthermore, the distribu- 
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tions of surface brightnesses of the two subsamples are 
similar — a K-S test yields a probability of 84% for the 
null hypothesis that they stem from the same underly- 
ing distribution. Certainly, measurements of the very 
central surface brightness, which are possible only with 
high-resolution observations, would provide a more di- 
rect argument here. However, since both nucleated and 
non-nucleated dEs within a given magn itude range have 
simil ar surface brightness profiles (Bing geli &: CameronI 
Il991f) . their effective surface brightness and central sur- 
face brightness are closely correlated. 

(2) Only one single galaxy among our 39 bright 
dE(nN)s (2.5%) is bright enough to fall among 
ICote et all 's definition of a HSB dE. In contrast, 14 of 
our 121 bright dE(N)s (12%) would quahfy as HSB dE. 
Therefore, it appears highly unlikely that a significant 
number of dE(nN)s possess nuclei with similar relative 
brightnes ses as those of the d E(N)s that were not de- 
tected bv lBinggeli et all (|1985f). 

(3) None of the dEs in iCote et all 's own sample that 
were previously classified as non-nucleated, but have now 
been found to host a weak nucleus, actually are HSB dEs. 

(4) Since we are interested in the distributions of our 
subsampl es with res pect to density in the cluster, we 
translate ICote et all 's conjecture about the spatial dis- 
tribution of the dEs into one about the distribution with 
respect to density: if the different density distributions 
of bright dE(N)s and dE(nN)s (see above) would primar- 
ily be caused by a surface brightness selection effect, a 
significantly larger fraction of the high surface brightness 
objects should be located at lower densities as compared 
to the lower surface brightness objects. To test for this 
possible bias, we plot the mean r band surface bright- 
ness within the half-light aperture against local projected 
density for the combined sample of bright dE(N)s and 
dE(nN)s (Fig. [5]). No correlation is seen, ruling out that 
such a bias is present in our data. 

We point out that it might of course still be the case 
that most of the dE(nN)s host weak nuclei that are be- 
low the VCC detection limit, as discussed in Sect. 12.21 
However, what is at stake here is the question whether a 
significant number of dE(nN)s should already have been 
classified as dE(N)s by the VCC, and whether this could 
account for the population differences that we find. The 
above arguments clearly rule out such a bias. We can 
thus conclude that the bright dE(N)s and dE(nN)s are 
indeed distinct dE subpopulations that differ in their 
clustering properties (Figs. [S] and [7]), as well as in their 
shapes (Figs. |3] and |4]). 

6. COLOR ANALYSIS 

Since our morphological subdivision of the dEs into 
several subpopulations is now established, the next step 
would obviously be to compare their stellar population 
properties. Given that the SDSS imaged every galaxy 
in five bands, it should be able to provide some insight 
into their stellar content, even though it is basically im- 
possible to disentangle ages and metallicities with optical 
broadband photometry alone. However, the issue is com- 
plicated by the fact that the u and z band images, which 
would be very important for an analysis of the stellar 
content, have a very low S/N (see Sect. [3]). It is there- 
fore important to perform a thorough study of the dE 
colors and color gradients that properly takes into ac- 



count measurement errors and the different S/N levels 
for objects of different magnitudes and surface bright- 
nesses. Such a study is beyond the scope of the present 
Paper and will be presented in a forthcoming Paper of 
this series. 

Nevertheless, in order to tackle the question about 
whether the dE subsamples differ in their color proper- 
ties, we present in Fig.[9]the inner u — g ("age-sensitive") 
versus i — z { "metallicity-sensitive" ) colors for the bright 
(mj. < 15.67) dEs, measured within an aperture of 
a — 0.5ahi,r- This approach guarantees relatively small 
errors (typical values are shown in the lower left corner 
of the Figure) that need not be taken into account indi- 
vidually. For each dE subsample, we indicate its median 
color values with the black symbols drawn with thick 
lines. 

However, a direct comparison of these values would be 
biased by the existence of a color-magnitude relation: if 
different subsamples had, on average, significantly differ- 
ent magnitudes, they would be offset in our color-color 
diagram even if they followed exactly the same color- 
magnitude relation. We therefore compute an approxi- 
mate correction for this effect: first, we perform a linear 
least squares fit to the color-magnitude relations (r ver- 
sus u — g and r versus i — z) of our full dE sample, clip- 
ping one time at 3(t and excluding the dE(bc)s because 
of their blue inner colors. We then derive the median 
r magnitude of each subsample, and use the linear fit 
to compute its expected color offset from the sample of 
dE(nN)s, which we choose as reference. The so obtained 
corrected median colors are shown in Fig. [5] as black sym- 
bols drawn with thin lines, and are connected with lines 
to their uncorrected values. 

The dE(bc)s exhibit, as expected, the bluest colors of 
all subsamples, basically by definition, since we focus 
here on the inner galaxy colors. While the corrected col- 
ors of the dE(di)s are similar to those of the dE(nN)s, the 
dE(N)s are, on average, redder in i — z and significantly 
redder in u — z. Given the very small color correction and 
large sample size of the dE(N)s, this can be considered a 
robust result. In the inset shown in Fig. [9l we compare 
the median values of the dE(nN)s and dE(N)s to two 
model tracks from stellar popu lation synthesis calcula- 
tions (|Bruzual fc Chariot! l2003f ). Both tracks represent 
stellar populations formed through a single burst of star 
formation that exponentially decays with time (r = 0.5 
Gyr), using Padova 2000 isochrones and a Chabrier IMF. 
The tracks are curves of constant metallicity (grey solid 
fine: Z = 0.008, grey dotted fine Z = 0.004); ages 
increase from bottom to top and are marked at 3, 6, 
10, and 14 Gyr (the latter mark is outside of the plot 
area for the Z — 0.008 track). Our measured values lie 
along of the Z — 0.008 track, illustrating that, within 
the framework of our simplified stellar population mod- 
els, the color difference between dE(nN)s and dE(N)s 
could be interpreted as a difference in age. According 
to this simple approach, the dE(N)s would be, on aver- 
age, a few Gyr older than the dE(nN)s. However, the 
measurements also fall roughly along a virtual line con- 
necting the 6 Gyr points of each model, showing that 
they might also be interpreted as a difference in metal- 
licity. While this color offset between the dE(N)s and 
dE (nN)s would be quali t atively consistent with the study 
bv iRakos fc Schomlpertl (j2004[ ). who find the dE(N)s in 
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the Coma and Fornax clusters to have older stellar pop- 
ulations than the dE(nN)s, reliable conclusions need to 
await a more comprehensive color study of our dEs. 

7. DISCUSSION 

7.1. Interrelations between subclasses 

The bright dE(nN)s and dE(di)s are both unrelaxed 
populations of relatively bright dEs shaped like thick 
disks. This also applies to the dE(bc)s, which could thus 
be candidates for being the direct progenitors of the for- 
mer: the presently blue centers of the dE(bc)s will evolve 
to typical dE colors within 1 Gyr or less after the cessa- 
tion of star formation (Paper II). Therefore, the bright 
dE(nN)s and dE(di)s could constitute those disk-shaped 
dEs where central star formation has already ceased. To 
test this hypothesis, we make the following considera- 
tions. There are 39 bright dE(nN)s, as well as 30 dE(di)s 
with mr < 15.67 mag, 7 of which are non-nucleated. This 
adds up to 69 "non-star-forming, disk-shaped dEs", 23 
(33%) of which are nucleated. Among the dE(bc)s there 
are 15 galaxies with rrip < 15.67 mag, 6 (40%) of which 
are nucleated. Thus, the fraction of nucleated galax- 
ies would be compatible with our hypothesis within the 
errors, with the caveat that n uclei might still form in 
the centers of some dE(bc)s I'see lOh fc Li^[200?il and Pa- 
per II), which would raise the nucleated fraction of the 
dE(bc)s. 

Now, 43% of the non-star- forming, disk-shaped dEs are 
dE(di)s, i.e., show disk features (not only an overall disk 
shape). If the dE(bc)s would contain the same fraction of 
galaxies that display disk features, we would expect 6.5 
such objects among the 15 dE(bc)s, with a standard devi- 
ation of 1.9. The observed number of 4 lies within 1.3 cr of 
the expected value and could thus still be reconciled with 
the above picture. However, since not only the dE(di)s, 
but also the bright dE(nN)s are disk-shaped, why do the 
latter not display disk features like the dE(di)s? This 
could either indicate a correlation between the presence 
of a significantly bright nucleus and the presence of disk 
substructure, or it could imply that there is more than 
one formation path towards disk-shaped dEs. 

7.2. Formation mechanisms 

If dEs originated from galaxies that fell into the clus- 
ter, how long a go co uld this infall have taken place? 
iConselice et al.l ()2001[ ) derived a two-body relaxation 
time for the Virgo dEs of much more than a Hubble time. 
Even violent relaxation, which could apply for the case 
of infalling or merging groups, would take at least a few 
cros sing times tcr, with tcr ~ 1.7 Gyr for the Virgo clus- 
ter (jBoselli fc GavazzH l2006f ) . Therefore, the majority 
of dE(N)s or their progenitors should have experienced 
infall in the earliest phases of the Virgo cluster (w hich 
is a rather yourig stru cture, see lBinggeli et al.lll987l and 
lArnaboldi et al.|[2004[ ). or they could have formed in dark 
matter halos along with the cluster itself. All other dE 
subclasses are largely unrelaxed populations, implying 
that they have formed later than the dE(N)s, probably 
from (continuous) infall of progenitor galaxies. Our color 
analysis in Sect. [6] would support this view, since it finds 
that the inner colors of the dE(N)s can be interpreted 
with an older stellar population than the dE(nN)s. This 
would be expected if one assumes that the progenitor 



galaxies had been forming stars until their infall into 
the cluster, resulting in a younger stellar population on 
average in the case of a later infall (neglecting possible 
metallicity differences). However, as stressed in Sect. |6l 
robust conclusions need to await a more detailed multi- 
color study of our dEs. 

The galaxy harassment scenario (|Moore et al.l Il996f) 
describes the structural transformation of a late-type spi- 
ral into a spheroidal system through strong tidal interac- 
tions with massive cluster galaxies. A thick stellar disk 
may survive and form a bar and spiral features that can 
be retained for s ome time, depending on t he tidal heat- 
ing of the galaxy ([Mastro2ietroIeF^ni2Q03). Harassment 
could thus form disk-shaped dEs, and dE(di)s in partic- 
ular. Moreover, it predicts gas to be funne led to th e cen- 
ter and form a density excess there (Mo ore et al.l[l998) . 
which would be well suited to explain the central star 
formation in the dE(bc)s. Therefore, it appears possible 
that harassment could form disk-shaped dEs that first 
appear as dE(bc)s and then passively evolve into dE(di)s 
and bright dE(nN)s as their star formation ceases (see 
Sect. 17. l() . It might also provide a way to form the fainter 
non-nucleated dEs, assuming that the tidal forces have a 
stronger effect on the shape of less massive galaxies, re- 
sulting in rounder objects on average. However, in order 
to explain all these subclasses by a single process, one 
would need to invoke a correlation between the presence 
of a nucleus and of disk fea tures, as discussed iii Sect. 17. II 

Ram-pressure stripping (jGunn &: Gottl[T"972[ ) of dwarf 
irregulars (dirrs) could be responsible for the fact that 
the disk-shaped bright dE(nN)s do not show disk fea- 
tures like the dE(di)s: dIrrs typically have no nucleus, 
and ram-pressure stripping exerts much less perturb- 
ing forces than a violent process like harassment, thus 
probably not triggering the formation of bars or spiral 
arms. Commonly discussed problems with this s cenario 
are the metallicity offset between dEs and dIrrs ( ThuanI 
Il985t [Richer et al.]|1998l: iGrebel et al][2003 ) and the too 
strong fading of dIr rs after cessation of star formation 
(jBothun et al.lfl986h . Also, the flattening distribution of 
Virgo cluster dIrrs - with intrinsic (primary) axial r a- 
tios > 0.5 for most galaxies (Binggeli fc Popesciilll995l ) - 
is not quite like that of our bright dE(nN)s. On the 
other hand, significant mass loss due to stripped gas 
might affect the stellar configuration of the galaxies and 
could thus possibly account for the difference. Moreover, 
the flattening distribution of the dIr rs is similar to that 
of th e faint dE(nN)s (cf. Fig. 9 of iBinggeli fc Popescul 
Il995f ). suggesting that these - and possibly not the bright 
dE(nN)s - might be stripped dIrrs. 

Tidally induced star formation of dIrrs might be able to 
overcome the problems of the ram-pressure stripping sce- 
nario: the initially lower metallicity and surface bright- 
ness of a dirr are increased by sev eral bursts of star 
formation (jPavies fc Phillippi Il988l ). during which the 
galaxy appears as blue compact dwarf (BCD). After the 
last BCD phase it fades to become a dE, thereby provid- 
ing an explanation how BCDs could be d E progenitors, 
which has frequently been discussed (e.g . , iBothun et al.l 
figsllPaoaderos et al.lfT99llGrebellfT997t Paper II). The 
last star formation burst might occur in the central re- 
gion, consistent with the appearance of the dE(bc)s. 

In addition to the number of possible formation sce- 
narios, the role of the nuclei provides another unknown 
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element. If dE(N)s and dE(nN)s would actually form a 
continuum of dEs with respect to re l ative n ucleus bright- 
ness as suggested by iGrant et al.l ()2005D . their signif- 
icantly different population properties could be inter- 
preted with a correlation between relative nucleus bright- 
ness and host galaxy evolution. Such a correlation could, 
for example, be provided by nucleus formation through 
coalescence of globular clusters (GCs): the infall and 
merging of several GCs - resulting in a rat her bright 
nucle us like in a dE(N) - takes many Gyr (|0h &: LinI 
|2000|) . consistent with the dE(N)s being in place since 
long. The dE(nN)s, on the other hand, were probably 
formed more recently, leaving time for only one or two 
GCs, or none at all, to sink to the center. 

7.3. Remarks on previous work 

Results similar to ours were derived by 

guson fc Sandagd (|1989( ). who also subdivided 
Virgo and Fornax cluster dEs with respect to magnitude 
and the presence or absence of a nucleus. In accordance 
with our results, they found that the dE(N)s are 
centrally clustered like E and SO galaxies, while the 
bright dE(nN)s are distributed like spiral and irregular 
galaxies. They also found the axial ratios of the bright 
dE(nN)s to be flatter than those of the dE(N)s. 

However, despite these similar results, their magni- 
tude selection of "bright" and "faint" subsamples is ac- 
tually quite different from ours. We initially selected 
only dEs with me < 18.0 mag (the completeness limit 
of the VCC), yielding a sample range of about 4.5 mag 
in B, and then subdivided our full sample at its median 
r mag nitude. In contrast to that, iFerguson &: Sandagd 
(|1989f ) included VCC galaxies with tub < 17.5 in their 
bright subsample, which therefore still spans a range of 4 
mag. Their faint subsample contains VCC galaxies with 
tob > 18.4, which are not included in our study and al- 
ready lie within the lu minosity regime of Local Group 
dwarf spheroidals (e.g.. iGrebel et all 120031 ). Therefore, 
their and our study can be considered complementary to 
some extent, in the sense that we probe different lumi- 
nosity regimes with our respective subsample definitions. 

8. CONCLUSIONS 

We have presented a quantitative analysis of the in- 
trinsic shapes and spatial distributions of various sub- 
samples of Virgo cluster early- type dwarfs (dEs) : bright 
and faint (non-)nucleated dEs (dE(N)s and dE(nN)s), 
dEs with disk features (dE(di)s), and dEs with blue cen- 
ters dE(bc)s). The dE(bc)s, dE(di)s, and bright dE(nN)s 
are shaped like thick disks and show basically no cen- 
tral clustering, indicating that they are unrelaxed pop- 
ulations that probably formed from infalling progenitor 
galaxies. As opposed to that, the dE(N)s (both bright 
and faint) are a fairly relaxed population of spheroidal 
galaxies, though an oblate intrinsic shape is favored for 
them as well. The faint dE(nN)s appear to be somewhat 
intermediate: their shapes are similar to the dE(N)s, but 
they form a largely unrelaxed population as derived from 
their clustering properties. Taken together, these results 
define a morpholog y-density relatio n within the dE class. 

Given that iFerguson fc Sandagd (|T98a ) derived similar 



results for both Virgo and Fornax cluster galaxies, it is 
also clear that this zoo of different dE subclasses is not 
only specific to the Virgo cluster. Similarly, a significant 
number of Coma cluster dEs show a two-c omponent pro- 
file a nd are flatte r than the normal d Es ( Agucrri et al.l 
120051 ). Moreover. iRakos fc SchombertI (|2004 ) found the 
dE(N)s in Coma and Fornax to have older stellar popu- 
lations than the dE(nN)s, consistent with our color anal- 
ysis of the Virgo cluster dEs. Thus, although the relative 
proportions of the dE subclasses might vary between the 
dynamically different Virgo, Coma, and Fornax clusters, 
the dE variety itself is probably similar in any galaxy 
cluster in the present epoch. We thus consider it im- 
portant that future studies of dEs do not intermingle the 
different subclasses, but instead compare their properties 
with each other, e.g., their stellar content or kinemati- 
cal structure. This will eventually lead to pinning down 
the actual significance of the various suggested forma- 
tion paths, thereby unveiling an important part of galaxy 
cluster formation and evolution. 
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Fig. 1. — Nucleus detection limits. Shown are central surface brightnesses and nucleus magnitudes in g, both from 'Cote et al.' f20 03), for 
34 of t he 36 VCC dEs of their sample (VCC 1512 has no nucleus, and for VCC 1743 no values could be derived). Objects where Cot e et al.l 
ll200g) identified "new" nuclei, i.e., that are listed as non-nucleated or only possibly nucleated in the VCC, are shown as filled squares. 
Objects that might have an offset nucleus according to Cote et al. (2006) and that are listed as non-nucleated in the VCC are shown as 
open circles. Objects classified as nucleated in both iCote et al.. 1,2006. ) and the VCC (including class "N:") are represented by the filled 
circles. The dotted line follows equal values of central surface brightness (in mag arcsec"^) and nucleus magnitude (in mag). 



Eorly-type dwarf ("dE", B^18) 



normal 

(83-88%) 



disk feature blue center 
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Fig. 2. — The zoo of early-type dwarfs. Our subdivision is shown as tree scheme, along with a 3-D illustration of a typical intrinsic 
shape, using the median axial ratio from the intrinsic distributions shown in Fig. [3] (see Sect. [5T2t . We also provide a statement about the 
inferred dynamical status of each subpopulation, as indicated by their clustering properties: "relaxed" and "unrelaxed" (see Sect. l5.3)l . The 
subdivision of dE(N)s into faint and bright samples is shown in grey color to reflect the fact that we find them to be not different in their 
properties (see Sect.[5ll. The branch of the dE(bc)s is shown as dotted line only, since these are not a morphological subclass (see text). The 
percentage of each subsample among our 413 Virgo cluster dEs is given in parentheses. For the dE(bc)s and dE(di)s the percentage ranges 
include corrections for the estimated number of objects missed by our detection techniques (see Papers I and II). Three sample images are 
shown for the normal dEs. For the dE(di)s, we show one sample image along with its unsharp mask revealing the spiral substructure. For 
the dE(bc)s, one sample image is shown along with its g — i color map revealing the blue center (dark=blue). See text for further details. 
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Note. — Classifi cation of a dE as nucleated or non-nucleated is provided by the VCC 
l IBinggeli et al.lll985D . Of the dEs that do not display disk substructure or a blue center, 
those with a small uncertainty on the presence of a nucleus ("N:") were included in the 
nucleated subclass, while those with a larger uncertainty ("N?", "Npec") were not assigned 
to any subclass (entry " — " as subclass), but were excluded from all comparisons of dE 
subclasses. Objects VCC 0218, 0308, 1684, and 1779 are dE(bc)s with disk features. 
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Fig. 3. — Deduction of intrinsic shapes. In the leftmost column, the number of objects in each of our dE subsamples are given. In the 
second column, we show the distribution of projected axial ratios for each of our dE subsamples. The data are shown as running histogram 
with a bin width of 0.15, corresponding to one fifth of the range in axial ratio covered by our galaxies. Each curve is sampled in steps of 
0.04 (one quarter of the bin width) and normalized to an area of 1. The observed distributions are approximated by analytic functions (see 
text), shown in the third column. Prom these, we derive intrinsic axial ratio distributions (fourth column), adopting purely oblate (grey 
solid line) and purely prolate (black dash-dotted line) shapes. For these distributions, we show in the fifth column 3-D illustrations of the 
galaxy shapes for the oblate (upper) and prolate (lower) case. For each distribution, we show the shape when using the 25th percentile 
axial ratio (left) and the 75th percentile axial ratio (right). In the rightmost column, we show for each dE subsample the surface brightness 
test (see text); we plot the surface brightness offset from the mean relation of r band surface brightness and magnitude against axial 
ratio. Surface brightness is measured within a = 2ahl,r> since axial ratio is measured at the same semimajor axis. The mean relation of 
surface brightness and magnitude is obtained through a linear least squares fit with one 3iT-clipping. The arrows pointing from the surface 
brightness test diagrams towards the shape illustrations reflect whether the test implies oblate or prolate intrinsic shapes; see text for more 
details. 
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Fig. 4. — K-S test results for the comparison of the axial ratio distributions of our dE subsamples. For each pair of distributions we give 
the probability in percent for the null hypothesis that the two distributions stem from the same underlying distribution function. 



Virgo early-type dwarfs. III. Subpopulations 



15 



dE(N) 
bright 

121 obj. 

(83 obj.) 

dE(N) 
faint 

89 obj. 

(32 obj.) 

dE(nN) 
faint 
110 obj. 
(13 obj.) 

dE(nN) 
bright 
39 obj. 
(19 obj.) 

dE(di) 
33 obj. 
(28 obj.) 

dE(bc) 
17 obj. 
(15 obj.) 



15 
10 

15 
10 

15 

10 



o 

o 15 

CN 

^2-10 
o 

UJ 

o 



15 
10 

15 
10 



"T": — r 



• >. • 



-HH- 



> - 



L V-i -J 



-I-H- 



190 185 
RA (J2000) 




2000 
Velocity / km s"' 



Fig. 5. — Spatial and velocity distribution. In the left column, the number of objects in each of our dE subsainples are given. Numbers 
in parentheses apply to those dEs for which heliocentric velocities are available. The middle column shows the projected distribution of the 
dE subsamples (black dots) within the cluster. All Virgo cluster member galaxies are shown as small grey dots. The right column shows 
the velocity distributions of the dE subsamples. The data are shown as running histogram with a bin width of 384km/s, corresponding 
to the semi-interquartile range of the total 193 velocities. Each curve is sampled in steps of 96km/s (one quarter of the bin width) and 
normalized to an area of 1. 
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Fig. 6. — Morphology vs. de nsity. Cumulat i ve dis tribution of local projected densities of our dE subsamples, and of Hubble types (inset). 
Following Dressier (1983) and lBinegeli et al.l ||T987^ . we define a circular area around each galaxy that includes its ten nearest neighbours 
(independent of galaxy type), yielding a projected density (number of galaxies per square degree). 
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Fig. 7. — K-S test results for the comparison of the distributions of our dE subsamples with respect to local projected density. For each 
pair of distributions we give the probability in percent for the null hypothesis that the two distributions stem from the same underlying 
distribution function. 
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Fig. 8. — Surface brightness vs. density. The mean surface brightness in r within the half-light aperture is compared to loc al projected 
densit y for the combined sample of bright dE(N)s and dE(nN)s, in order to test for a possible classification bias as conjectured bv lCote et al.l 
II2006I) (see text). No correlation is seen, ruling out such a bias. 




0.1 0.2 

Fig. 9. — Distribution in color space. Shown are the inner u — g versus i — z colors, measured within a = O.Sajji ^, for all dEs brighter than 
the median r brightness of our full sample, rrir = 15.67 mag, divided into the respective subsamples. Individual measurements are shown 
with small grey symbols (see the legend above the diagram). The median value of each subsample is shown as black symbol using thick 
lines. Black symbols drawn with thin lines represent the median values corrected for the effect of the color-magnitude-relation (see text); 
lines connect them to the corresponding uncorrected values. The correction is chosen to be zero for the dE(nN)s. The inset shows again 
the median values (corrected and uncorrected) of the dE(N)s and dE(nN)s, along with two model tracks from stellar population synthesis 
calculations (Bruzual & Chariot 2003). Both tracks represent stellar populations formed through a single burst of star formation that 
exponentially decays with time (r = 0.5 Gyr), using Padova 2000 isochrones and a Chabrier IMF. The grey solid line is for a metallicity 
Z = 0.008; age steps are marked by the grey diamonds at 3 Gyr, 6 Gyr, and 10 Gyr. The grey dotted line is for Z = 0.004, with age steps 
marked by the grey circles at 3 Gyr, 6 Gyr, 10 Gyr, and 14 Gyr; the latter is also the end of the track. 



